Impurities block the a to to martensitic transformation in titanium 
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The onset and kinetics of martensitic transformations are controlled by impurities trapped during 
the transformation. For the a-^ui transformation in Ti, ab initio methods yield the changes in both 
the relative stability of and energy barrier between the phases. Using the recently discovered trans- 
formation pathway, we study interstitial O, N, C; substitutional Al and V; and Ti interstitials and 
vacancies. The resulting microscopic picture explains the observations, specifically the suppression 
of the transformation in A- 70 and Ti-6A1-4V titanium alloys. 

PACS numbers: 



Impurities control the onset and kinetics of martensitic 
phase transformations — diffusionless structural transfor- 
mations proceeding near the speed of sound [1|. The 
abundant use of martensitic materials in engineering 
technologies is not underpinned by a microscopic under- 
standing of the effects of impurities during the trans- 
formation. The a— >lu martensitic transformation in Ti 
lowers its toughness and ductility, affecting its use in the 
aerospace industry. The 2 to 9 GPa range in the pub- 
lished transition pressures of nominally pure Ti 0, Q is 
believed to be caused by sample purity 0, 0] . System- 
atic studies show a large sensitivity to small amounts of 
impurities 0, • Two commercial Ti alloys — A- 70 and 
Ti-6A1-4V — show no transformation up to 35 GPa 0. 

Impurities pose two theoretical challenges: The effect 
on the relative stability of the two phases and the energy 
barrier of the transformation. Understanding both of 
these requires the atomistic pathway for a— hj to model 
the impurity motion during the transformation. The near 
sonic speed of the transition prohibits diffusion during 
the transformation, thus requiring impurities to move 
collectively with their local environment. Recently Trin- 
kle et al. identified from all possible transformation path- 
ways the a— >lu mechanism with the lowest energy barrier 
of 9 meV/atom. 

This paper exploits this fast transformation that traps 
the impurities in their known path to systematically 
study a range of impurities, including those in the com- 
mercial Ti alloys A- 70 and Ti-6A1-4V. (1) We study inter- 
stitial O, N, C; substitutional Al and V; Ti interstitials 
and vacancies. These impurities and defects affect the 
transformation by shifting the relative stability of and 
energy barrier between the a and uj phases. (2) Intersti- 
tial impurity effects are governed primarily by their size, 
while substitutional impurities affect the transformation 
by changing the d-electron concentration. (3) The most 
important impurities, O in A-70 Ti and Al in Ti-6A1-4V, 
more than double the transition barrier and increase the 
energy of u> relative to a, explaining the observed sup- 
pression of the a^uj transition. 

Method - The ab initio calculations are performed with 



VASP [8J, |9j, a density functional code using a plane- 
wave basis and ultrasoft Vanderbilt type pseudopoten- 
tials 0, . The generalized gradient approximation of 
Perdew and Wang is used [12]. A plane- wave kinetic- 
energy cutoff of 400 eV ensures energy convergence to 
1 meV/atom. The k-point meshes for the different struc- 
tures are chosen to guarantee the same accuracy. To 
avoid unphysical short distances between the atoms, it 
is necessary to include the Ti 3p states as valence states 
in addition to the usual 4s and 3c? states. For the other 
elements the pseudopotentials describe the core states as 
follows: V [Ar], Al [Ne], O, N and C [He]. 

Impurity locations and formation energies are deter- 
mined by relaxations using 4x4x3 and 3x3x4 su- 
percells for the a and ui phases with 96 and 108 atoms, 
respectively. This results in a 1% impurity concentra- 
tion. A 2 x 2 x 2 k-point mesh is used. The atom 
positions are relaxed until the total energy changes by 
less than 1 meV, corresponding to atomic-level forces 
F max < 20meV/A. Altogether this provides an energy 
accuracy of 0.5 meV/atom. 

NEB calculations with variable cell shape at constant 
pressure yield the energy barriers of the martensitic 
transformation [l3 |. The pathway is represented by 18 in- 
termediate images of a 48 Ti atom supercell constructed 
using y/2 x 2^2 x 2 cells of the TAO-1 mechanism with 
a 4 x 4 x 4 k-point mesh. This results in an impurity 
concentration of 2%. Comparing the formation energies 
for the 48 atom supercell with the larger cells of 96 and 
108 atoms provides an estimate of the finite size error 
of 0.05 eV. The images are relaxed until the energy of 
each image changes by less than 1 meV, corresponding 



to atomic-level forces F„ 



< 50meV/A and stresses 
Cmax < 20MPa. Altogether this yields an energy bar- 
rier accuracy of 1 meV/atom, the same as the finite size 
error. To compare the effects of different defects, all en- 
ergy barriers are reported relative to the a structure for 
that defect, and in units of meV per 48 Ti atoms. 

Interstitial impurities - Figure ^ shows the location of 
interstitial impurities and the effect of the transformation 
on these sites. The interstitial impurity sites were found 
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FIG. 1: Octahedral (o), tetrahedral (t), and hexahedral (h) 
sites for interstitial impurities as well as A and B sites for 
substitutional impurities in the a and u> phases. The a and uj 
phases each contain one unique o, t, and h site. The a ta t site 
relaxes to the nearby «hox site for all three impurities (O, N, 
C). The hexahedral site is a distorted double-tetrahedral site 
with five neighbors. The arrows indicate the transformation 
of the impurity and lattice sites in the TAO-1 mechanism, 
and the relative proportions. 



by placing O, N, and C impurities on all interstitial sites 
in a and uj and relaxing the atomic positions. Each phase 
contains one unique octahedral, tetrahedral, and hexahe- 
dral site. However, the tetrahedral site in a is unstable 
for O, N, and C and relaxes to the nearby hexahedral site 
in the basal plane. The "hexahedral" site is a distorted 
double-tetrahedral site with five neighbors, formed by a 
triangle of the basal plane and the two atoms right above 
and below the center of the triangle. The instability of 
the tetrahedral site is unexpected from crystallographic 
considerations [l4| and is due to atomic relaxations. In 
uj the tetrahedral and hexahedral sites are both stable, 
and they are close in energy and location. 

Table |Ja) shows the formation energies of the inter- 
stitial defects in a and uj. The energy is measured rela- 
tive to molecular O2, N 2 and graphite, respectively. For 
both a and uj the octahedral site is more stable than the 
tetrahedral or hexahedral sites. This is simply related to 
the size of the interstitial sites: larger sites have lower 
formation energies. Comparison between the a and uj in- 
terstitial impurities shows that the octahedral formation 
energy is slightly lower for a; this difference shifts the 
relative energy of the two crystal structures, decreasing 
the stability of uj over a. 



Figures Ufa) and (b) show octahedral sites in a — 
favored by O, N, and C — transforming into two possible 
sites in uj and increase the energy barrier for the transfor- 
mation. The TAO-1 pathway breaks the hexagonal sym- 
metry and transforms a third of the octahedral sites into 
the uj octahedral site and two thirds into the uj hexahedral 
site. The presence of impurities increases the energy bar- 
rier regardless of the final site of the impurities; however, 
this increase is much larger for the octahedral to hexahe- 
dral pathway than for the octahedral to octahedral one. 
For all three impurities, despite their different chemistry, 
their similar size results in nearly identical barriers. 



TABLE I: Formation energies and location of interstitial 
and substitutional impurities as well as vacancies and self- 
interstitials in the a (PG^/mmc) and uj (P6/mmm) phases. 
The formation energies Ef for O, N, and C are measured rel- 
ative to molecular O2, N2, and graphite, respectively, and for 
Al and V relative to their fee and bec phase, respectively. For 
the vacancies and self-interstitials the formation energies are 
measured relative to the corresponding a or uj phases. All cal- 
culations are performed at 1 at.% defect concentration. For 
each site, the average nearest neighbor distance after relax- 
ation, Rnn, is similar for all impurities; hence, we show the 
range of values. The coordination number Z at each site is 
also listed. 



Site Wyckoff pos. a R NN [A] Z 



Ei [eV] 



(a) Interstitial impurities 
a oct 2(a) (0,0,0) 2.06-2.09 
2(d) a 



Oh 



'2 11' 

v 3 ' 3 ' a 



1.91-1.95 



Woct 3(f) (|,0,0) 
^hcx 6(k) (x,0, §) 
ujtct 6(m) (a;, 2x 



2.03-2.17 
1.92-2.18 
1.90-2.02 



O 
-6.12 
-4.93 

-6.06 
-4.47 
-4.38 



N C 

-5.10 -1.58 

-3.41 +0.49 

-4.99 -1.36 

-3.03 +0.61 

-2.83 +0.78 



(b) Substitutional impurities 

a 2(c) (§,§,£) 2.85-2.92 12 

uja 1(a) (0,0,0) 2.82-3.00 14 

lub 2(d) (1 2 ,|) 2.60-3.00 11 



Al V 

-0.88 +0.51 

-0.90 +0.60 

-0.47 +0.33 



(c) Vacancy 

a 2(c) 

uja 1(a) 

ujb 2(d) 



(1 2 I) 

V 3 ' 3 ' 4,1 
(0,0,0) 

(I 2 L) 

v 3 ' 3 ' 2 1 



2.85-2.94 12 
2.83-3.02 14 
2.12-3.03 11 



+2.03 
+2.92 
+1.57 



(d) Self-interstitial 



Q?oct 

Odumb 

UJacX, 

Whex 

UJtcX 



2(a) 
2(d) 
3(f) 
6(k) 
6(m) 



(0,0,0) 

(2 I I) 

\ 3 ' 3 ' 4> 
(|,0,0) 

(SB, 0,3) 

(x,2x,\) 



2.47 

2.05 
2.33-2.60 
2.21-2.55 
2.15-2.51 



+2.58 
+2.87 
+3.76 
+3.49 
+3.50 



"See Ref. [jj 

Substitutional impurities - Table Hlb) shows the for- 
mation energies of the substitutional Al and V for a and 
the A and B sites in uj. These energies are defined rela- 
tive to fee Al and bec V. Comparison between the a and 
uj substitutional impurities shows that, for both Al and 
V, the favored uj site is lower than the a site. However, 
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FIG. 2: Energy barrier for the TAO-1 mappings of intersti- 
tial and substitutional impurities, as well as vacancies and 
interstitials in Ti relative to the a phase. The defect con- 
centration is 2%. The endpoint energies match the formation 
energies of Table [I] for a defect concentration of 1 at.% within 
1 meV/atom, providing an accuracy estimate for the barrier. 
The interstitial O, N, and C impurities occupy the octahedral 
sites in a and transform either (a) into the octahedral site in 
uj or (b) into the hexahedral site, with a 1 : 2 relative ratio. 
The substitutional Al and V impurities transform in uj either 
(c) to the A site or (d) to the B site, with a 1 : 2 relative ratio. 
Figure (e) shows the barrier for vacancies and (f) for the Ti- 
interstitial defects. The vacancy transforms from the a sites 
to the A and B sites in uj with a ratio of 1 : 2. The interstitial 
transforms from the octahedral site in a to the octahedral 
and hexahedral sites in uj, also with a ratio of 1 : 2. 

during the transformation, a random mixture of uj oc- 
cupations will be produced in a 1 : 2 ratio as shown in 
Figure^] In that case, the formation energy combination 
wa/3 + 2wb/3 is lower than a for V, and higher for Al. 
This is expected to increase the stability of uj relative to 
a with V, and decrease the stability of uj with Al. 

Figure|2c) and (d) show how the a substitutional sites 
for Al and V transform to the two possible uj sites and 
effect the transformation barrier. The TAO-1 pathway 
lowers the symmetry and transforms a third of the a 
sites into uja sites and two thirds into wb sites. The 
presence of impurities changes the energy barrier in a 
similar fashion as the formation energies. Since the shift 
in the barrier matches the shift in the relative energies of 
a and uj, we conclude it is due to the d-electron concen- 
tration. The ci-electron effect is also seen in the alloying 
behavior: Transition metals with large ci-electron con- 



centrations such as V, Mo, Fe, Ta stabilize the uj phase, 
while simple metals and early transition metals such as 
Al favor a 0. 

Titanium vacancies and interstitials - Table QJc) and 
(d) show the formation energies of the vacancies and Ti- 
interstitial defects in a and uj, respectively. The reference 
energy is the appropriate Ti crystal phase. For the vacan- 
cies, the ujb site — which relaxes the 6-fold B-sublattice 
rings into 5-fold rings — is most forgiving. However, dur- 
ing the transformation, a random mixture of uj occupa- 
tions will be produced in a 1 : 2 ratio; in that case, the 
formation energy combination wa/3 + 2ojb/3 is nearly 
identical to the a vacancy. The interstitials are energeti- 
cally favorable, as is expected for a closed-packed metal. 
While the a-octahedral site has the lowest formation en- 
ergy, there is an energetically close defect formed when 
the tetrahedral interstitial relaxes and forms a [0001] 
dumbbell configuration. All the uj interstitial sites are 
nearly degenerate in energy. The increased formation 
energy for uj interstitials over a interstitials is expected 
to shift the relative energy of the two crystal structures, 
decreasing the stability of uj over a. 

Figure [^e) and (f) show how vacancies and Ti- 
interstitials transform to two possible sites in uj and in- 
crease the energy barrier for the transformation. The 
TAO-1 pathway acts on vacancies and Ti-interstitials as 
for substitutional and interstitial impurities, respectively. 
The presence of vacancies increases the barrier in both 
cases, though the effect is much more pronounced when 
transforming to the uja vacancy. The presence of inter- 
stitials increases the energy barrier regardless of the final 
site of the interstitial; unlike the impurities, there is not 
much difference between the octahedral and hexahedral 
pathways. Due to their inherently low concentration the 
self defects will not have any significant effect on the 
transformation. 

Discussion - The a-uj phase transition is sensitive to 
small energy changes. At zero pressure and tempera- 
ture, the uj phase is about 4 meV/atom lower in en- 
ergy than a 1 161 and the energy barrier from a to uj 
is 9 me V/ atom |l7|. It takes 10 GPa to overcome the 
energy barrier in shock experiments or 280 K to over- 
come the energy difference [l^. Hence, we expect a sig- 
nificant change in transition temperature and pressure 
due to the energy shift caused by impurities. 

Table ITT1 shows the effect of impurities in A- 70 Ti and 
Ti-6A1-4V, especially the significant effect of O and Al. 
Overall, we find the presence of impurities in both al- 
loys suppresses the a— >uj martensitic phase transition by 
increasing both the uj energy relative to a and the en- 
ergy barrier. The increased energy barrier is expected to 
reduce nucleation and to slow down the transformation. 
The change of the a-uj energy difference shifts the a-uj 
phase boundary 0] downward by several hundred de- 
grees and thus at room temperature increases the trans- 
formation pressure dramatically. 
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TABLE II: Change of the relative energy between a and u> 
and of the energy barrier by impurities in the A- 70 and the 
Ti-6A1-4V alloys. The impurity concentrations are taken from 
Ref-B 



Alloy 


Tin ni i n tv 


A E I'm e VI 


AR, fmeVl 


A- 70: 


O (1.10 at.%) 


+12 


+10 




N (0.08 at.%) 


+ 1 


+1 




C (0.07 at.%) 


+ 1 


+1 




Total 


+14 


+ 12 


Ti-6A1-4V: 


Al (10.7 at.%) 


+29 


+31 




V (3.8 at.%) 


-3 


-3 




O (0.5 at.%) 


+6 


+5 




Total 


+33 


+33 



The other interstitial impurities, N and C, have an ef- 
fect similar to O. Despite their different chemistry, all 
increase the energy barrier and the relative energy by 
nearly the same amount. This indicates the primary ef- 
fect is steric, and hence, we expect other interstititals 
of comparable size to have similar effects. In the A-70 
Ti and Ti-6A1-4V alloys the effect of N and C is smaller 
proportional to their lower concentration. 

In contrast, the substitutional Al and V impurities 
have opposite effects due to the change in d-electron con- 
centration. As expected from the alloying behavior of 
Ti 0, V decreases the energy of uj and Al increases it. 
We observe the same behavior for the energy barriers. 
The Al impurity reduces the d-electron concentration by 
two, while V increases it by one; thus Al has a larger 
effect on the a— >uj transformation. Commercial Ti-6A1- 
4V alloy contains both; however, the concentration of Al 
(11 at.%) is three times higher than the concentration of 
V (4 at.%), explaining the observed suppression of the 
a— no transformation. 

Conclusion - We determined the energy and location 
of point defects in a and uj Ti and showed how they sup- 
press the martensitic transformation. Interstitial O, N, 
and C occupy the octahedral interstitial sites in a and 
uj and transform into the octahedral and hexahedral in- 
terstitial sites in uj. The interstitial impurities retard 
the transformation by increasing the transformation en- 
ergy barrier while the substitutional impurities Al and V 
influence the barrier by changing the (i-electron concen- 



tration. The effect of impurities on relative energies and 
energy barriers is central to understanding martensitic 
transformations. 
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